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Gamma-ray bursts (GRBs) and their afterglows have been proposed as an excellent probe 


to study the evolution of cosmic star formation, the reionization of the intergalactic medium, 
and the metal enrichment history of the universe since the prompt 7-ray emission of 
GRBs should be detectable^ out to distances 2 ^ > 10. Hitherto, the highest measured redshift 
for a GRB has been = 4.50^. Here we report the optical spectrum of the afterglow of GRB 
050904 obtained 3.4 days after the burst. The spectrum shows a clear continuum at the long 
wavelength end of the spectrum with a sharp cutoff at around 9000 A due to Lya absorption 
at a redshift of 6.3 with a damping wing. Little flux is present in the waveband shortward of 
the Lya break. A system of absorption lines of heavy elements at redshift 2 ^ = 6.295 ± 0.002 
were also detected, yielding a precise measurement of the largest known redshift of a GRB. 
Analysis of the Si II fine structure lines suggest a dense metal-enriched environment around 
the GRB progenitor, providing unique information on the properties of the gas in a galaxy 
when the universe was younger than one billion years. 

GRB 050904 was a long burst (duration T 90 = 225 s) detected by the Swift 7 -ray burst 
satellite on 4 September 2005, 01:51:44 UT^’^. Its position was immediately disseminated via the 
GRB Coordinates Network (GCN). While the optical observations at the Palomar 60" telescope 
carried out 3.5 hours after the trigger did not reveal a new source with upper limits of i? > 20.8 
mag and / > 19.7 mag^, a relatively bright near infrared source with J ~ 17.5 mag was detected 
three hours after the burst in the Swift X-ray telescope error circle^, which showed a temporal 
decay with an index of —1.20, fully consistent with being a typical GRB afterglow. Analysis of the 
near-infrared colors, combined with the non-detection in the optical bands, led to the suggestion 
that the burst originated at a high redshift^*’, 5.3 < 2 : < 9.0. A refined photometric redshift 
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2 : = 6.10lgj2 was reported based on ESO VET observations in the J, H, K, and I bands^^ 

We observed the field of GRB 050904 with the Eaint Object Camera And Spectrograph 
(EOCAS)^^ on the 8.2m Subaru Telescope atop Mauna Kea starting on the night of 6 September. 
In the z' band image (600 s exposure, mid-epoch 7 September, 8:04 UT), we detected the afterglow 
at 2 ;'(AB) = 23.71 ±0.14 mag, but we failed to detect it in the Ic band even with a longer exposure 
(900 s), which implied that the Eyman break should be present around A ~8500-9000 A. 

We then obtained a grism spectrum of the afterglow candidate, which exhibited a sharp cutoff 
at A ~9000 A with strong depletion of the continuum at shorter wavelengths, strikingly similar to 
the spectra of quasarsat 2 : > 6 except for the absence of a broad Eya emission line. The emission 
is very weak in the wavelength range shorter than 8900 A. In particular the flux is consistent 
with zero in the ranges 8500-8900 A and 7000-7500 A that extend shortward of the Eya and 
Ey/9 wavelengths for 2 : ~ 6.3. This is a clear signature of absorption by neutral hydrogen in the 
intergalactic medium at 2 : > 6, and marks the first detection of a Gunn-Peterson trough^"^ from an 
object other than high-;^ quasars'^. We also find weak emission features at ~7500-8300 A which 
are presumably leakage flux from the continuum emission that is also found in quasar spectra at 
similar redshifts. 

At the longer wavelength end of the spectrum is a flat continuum with a series of absorption 
lines, which we identify as S II, Si II, O I, and C II lines at a common redshift of z = 6.295 ± 
0.002. We believe that this is the redshift of the GRB host galaxy, since no other absorption line 
system was observed at a redshift consistent with that of the Eya break. This firm spectroscopic 
identification of the redshift breaks the previous record of GRB 000131 at 2;=4.50^. 
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From a closer examination of the absorption lines, we find that they are not saturated and can 
be used to estimate the column densities of the heavy elements as shown in Table 1. Using the stan¬ 
dard photospheric solar abundances'^ we obtain the metalicity of these elements as [C/H]=—2.4, 
[0/H]=—2.3, [Si/H]=—2.6, and [S/H]=—1.0, where logiVHi(cm“^) ~ 21.3 is assumed based on a 
DLA (damped Lyman a system) model for the Lya damping wing presented below. These values 
may not represent the typical abundances in the GRB host galaxy for several reasons. First, they 
are derived using only a single ionization state for each element. Depletion due to dust condensa¬ 
tion may modify the Si abundance in particular. And second, the spatial distribution of the heavy 
elements may be significantly different from that of hydrogen. It is possible that the heavy ele¬ 
ments are distributed only locally around the GRB source in a metal-enriched circumstellar shell, 
while the neutral hydrogen is distributed on a larger scale in or outside the host galaxy. 

Further analysis of the Si lines allows us to constrain the scale of the absorbing metals. Using 
the equivalent width ratio of the fine structure transition lines Si lU A 1264.7 and Si II A1260.4, the 
electron density can be constrained'^ as logne(cm“^) = 2.3 ± 0.7 for a reasonable temperature 
range of 10^ K < T < 10^ K. Combined with the column density and the abundance of Si derived 
above and assuming a hydrogen ionization fraction of 0.1, we obtain the physical depth of the 
absorbing system to be 0.4 pc with an uncertainty of a factor of ~ 10, reflecting the statistical 
errors and the possible temperature range. These fine structure lines have been found in GRB 
afterglow spectra whereas they have never been clearly detected in QSO-DLAs'^. This is 
consistent with a local origin for the absorption such as a metal-enriched molecular cloud in the 
star-forming region or a dense metal-enriched shell nebula swept-up by a progenitor wind prior to 
the GRB onset suggested for GRB 021004^'’^^ and GRB 030226^^. The column density of C II is 
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also consistent with the calculation for a carbon-rich Wolf-Rayet wind^"^. 


As shown in Fig 1, the Lya cutoff exhibits the signature of a damping wing redward of the 
Lya wavelength. This is the first detection of significant neutral hydrogen absorption at z > 6, 
and offers the first case for which one can explore the distribution of neutral hydrogen in the 
vicinity of a GRB, in the host galaxy, and/or in intergalactic space at such high redshifts. Such a 
study is difficult with high- 2 : quasars due to their enormous UV flux which ionizes the surrounding 
environment, and due to the presence of a strong Lya emission line. 

There are two possibilities for the nature of the absorber. It may be a DLA associated with 
the host galaxy, which has been observed in the afterglows of several GRBs at lower redshifts 
The other possibility is the neutral hydrogen in the intergalactic medium (IGM) left over from the 
pre-reionization era^. If the latter is the case, we can measure the neutral fraction of the IGM at 
2 : > 6 for the first time, giving crucially important information on the reionization history of the 
universe. We find that the wing shape can be reproduced either by DLA absorption (see inset of 
Fig. 1) or by the IGM. A comprehensive spectral fitting analysis is necessary to examine these 
possibilities, which is beyond the scope of this letter. 

With the detection of metal absorption lines, we have shown that GRBs are found in metal- 
enriched regions even at such an early phase of the Universe as 2 : > 6. It is, therefore, possible 
that one would detect the metal absorption lines even from GRBs originating from the metal-free 
first generation stars, as their environment may be self-polluted by pre-burst winds as suggested 
by the present observation. One can expect to obtain afterglow spectra with much higher quality 
for GRBs at even higher redshifts in the immediate future, considering that Swift is constantly 
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localizing faint GRBs, and that our spectrum was taken when the afterglow had faded by more 
than an order of magnitude since its first detection in the J band. Such future data will give us 
even better opportunities to probe the formation of stars and galaxies in the early universe. 
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Observed 

Equivalent 

Column 

Line 

Redshift 

Wavelength 

Width 

Density 

Identification 


(A) 

(A) 

log (cm~2) 



9041.0 ±0.8 

4.5 ± 1.0 

o o 

+ 

C IV A1548.2 

4.840±0.001 



— 

(N VA1238.8) 

(6.298±0.001) 

9055.9 ± 1.7 

1.7 ± 1.0 


C IV A1550.8 

4.840±0.001 



— 

(N V A1242.8) 

(6.287±0.001) 

9146.4 ± 1.8 

3.8 ± 1.1 


S II A1253.8 

6.295±0.001 

9188.7 ± 2.6 

6.1 ± 3.7 

16.201;:'' 

S II A1259.5 

6.295±0.002 

9195.9 ± 1.2 

8.3 ± 2.7 

14.291S;|; 

Si II A1260.4 

6.296±0.001 

9225.8 ± 1.8 

3.9 ± 1.1 

13.631S;;' 

Si 11* A1264.7 

6.295±0.001 

9499.1 ± 0.9 

10.3 ± 1.9 

15.85l»;i 

0 1 A1302.2 

6.295±0.001 

9737.2 ±1.1 

12.3 ± 2.4 

15.411»;» 

C II A1334.5 

6.296±0.001 


Table 1: List of the absorption lines detected in the spectrum of the optical afterglow of 
GRB 050904. The wavelengths and equivalent widths were derived by fitting a single 
Gaussian. The column densities of lines were estimated by the standard curve of growth 
analysis^®. The equivalent widths in the table are observed ones, i.e., not converted to 
the rest-frame. The quoted uncertainties are 1a statistical errors. Most of the absorption 
lines are consistent with being at a single redshift of z = 6.295 ± 0.002 within the statis¬ 
tical uncertainties. The absorption lines at A =9041.0 and 9055.9 A could be identified 
as N V AA1238.8,1242.8 if they are at a similar redshift as the other absorption lines. 
However, the derived redshifts of these two lines are significantly inconsistent with each 
other. Another possible identification is C IV AA1548.2, 1550.8 in an intervening system 
at z = 4.840, which we think is more likely. 


10 





Figure 1 Spectrum of the afterglow of GRB 050904 covering wavelengths 7000-10000 
A with a resolution R = A/AA ~ 1000 at 9000 A. It was taken with Subaru/FOCAS at mid¬ 
epoch 7 September, 12:05 LIT, 3.4 days after the burst for a total exposure of 4.0 hours. 
The abscissa is the observed wavelength converted to that in vacuum. The spectrum 
is smoothed to a resolution R = A/AA ~ 600 at 9000 A. The locations of the identified 
absorption lines (see Table ^) at z = 4.840 and 2 : = 6.295, as well as the wavelengths 
of Lya and Ly/3 at z = 6.295 are shown with vertical dotted lines. The one-sigma errors 
are shown with an offset of -1.0 at the bottom of the panel. In the inset, we show a 
model for the damping wing of Lya absorption with a neutral hydrogen column density 
logiVHi(cm“^) = 21.3 at a redshift = 6.3 by the solid line, overlaid on the observed spec¬ 
trum in the wavelenghth range of 8700-9500 A. The dotted line indicates the unabsorbed 
continuum model following a power-law (/^ oc function as typically observed for GRB 
afterglows. Note that only the red wing is relevant to the fit, since the emission blueward 
of Lya is absorbed by the IGM. 
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